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Abstract.  Using  the  LLNL  EBIT/RKTRAP  system  non-neutral  plasmas  of  highly  charged  ions 
were  produced  and  cooled  to  temperatures  around  one  Kelvin.  These  strongly  coupled  plasmas 
can  model  white  dwarf  astrophysical  plasmas  in  the  laboratory.  These  systems  may  also  have 
potential  application  to  quantum  computation.  The  experimental  results  from  the  last  operations 
of  the  trap  at  Livermore  are  discussed.  Molecular  dynamics  simulation  results  are  discussed  as 
a guide  to  past  and  future  experiments.  The  status  and  future  plans  for  RETRAP  at  LBNL's  88 
inch  Cyclotron  are  di.scussed. 


INTRODUCTION 

Strongly  coupled  plasmas  are  characterized  by  the  ratio  of  the  Coulomb  interaction 
between  particles  with  charge  q to  the  thermal  energy  of  the  particles  (k|jT)  being 
greater  than  one.  This  Coulomb  coupling  parameter: 


with  the  Wigner-Seitz  radius 

Where  n is  the  number  density  of  the  charged  particles.  If  T exceeds  -172  (for 
infinite  one  component  plasma)  the  plasma  crystallizes  and  is  refened  to  as  a 
“Coulomb  crystal”  or  for  finite  crystals  at  higher  Fs  “Coulomb  cluster”[l].  White 
Dwarf  stars  contain  plasmas  with  fully  stripped  ions  and  degenerate  electrons  that  can 
be  treated  as  strongly  coupled  non-neutral  plasmas.  Crystallization  and  other 
thermodynamic  properties  (chemical  layering,  etc.)  effect  the  luminosity  function, 
which  is  used  for  astrophysical  chronomctry[2,3]. 

Cold,  confined  and  relatively  numerous  highly  charged  ions  are  of  interest  for 
precision  fundamental  spectroscopic  measurements[4].  These  transitions,  coupled 
with  the  charge  scaling  of  motions  in  a trap,  provide  an  additional  “parameter  space” 
for  the  investigation  of  trap  based  quantum  coherence  and  computation  experiments. 
To  study  mixed  strongly  coupled  non-neutral  plasmas  and  to  prepare  for  precision 
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fine-structure  and  hyperfine-structure  laser  spectroscopy  on  ground  term  transitions  of 
highly-charged  ions,  a source  of  highly  charged  ions  (EBIT)  and  a cryogenic  Penning 
trap  for  capturing  and  cooling  these  highly  charged  ions  (RETRAP)  have  been 
developed  at  Lawrence  Livermore  National  Lab  [5,6,7], 

EXPERIMENT 

The  scheme  used  in  the  experiment  to  cool  highly  charged  ions  to  the  lowest 
realizable  temperature  via  sympathetic  cooling  with  laser  cooled  ions.  Be”^  ions  are 
used  for  both  their  low  mass  to  charge  ratio  (m/q=9)  and  the  reasonable  cooling  laser 
transition  (reachable  with  commercial  available  internally  doubled  ring  dye  laser).  Be^ 
ions  were  produced  in  a metal  vapor  vacuum  arc  (MeVVA)  source.  The  Be^  ions  were 
analyzed  by  a 90  degree  bending  magnet  (to  eliminate  other  charge  states  and  species 
produced  in  the  MEVVA)  and  focused  through  a series  of  ion  optics  into  a tube  at  a 
potential  near  the  MEVVA  extraction  potential  of  approximately  7.5  keV.  A high 
voltage  switch  lowers  the  potential  of  the  tube  to  near  ground,  “de-accelerating”  the 
Be"^  ions,  so  they  can  then  relatively  slowly  pass  into  a hyperbolic  trap  where  an 
auxiliary  “capture”  electrode  is  pulsed  to  confine  the  ions.  These  (-10^)  ions  were 
initially  cooled  to  near  room  temperature  by  tuning  (via  trap  ring-endcap  voltage)  their 
axial  oscillation  frequency  into  resonance  with  a cryogenic  parallel  tuned  circuit.  The 
heating  of  the  parallel  circuit  by  currents  induced  by  fluctuations  in  the  centroid  of  the 
confined  ion  motion,  together  with  ion-ion  collisional  coupling,  initially  cools  the  ions 
at  an  exponential  rate  with  a time  constant  of  several  minutes. 

Eventually,  the  cooling  is  dominated  by  a transverse  laser  beam  de-tuned  several 
GHz  below  the  cycling  resonance  transition  near  313  nm  [8].  With  decreased  de- 
tuning of  the  cooling  laser,  the  ions  then  rapidly  cooled  to  below  1 K,  as  evidenced  by 
a decrease  in  light  scattering,  and  confirmed  by  probing  the  Be  ions  with  a weak 
second  laser  tuned  through  a nearby  non-cycling  transition.  Using  slightly  altered  trap 
potentials  more  Be^  ions  can  be  ‘stacked’  into  the  trap,  with  a cooling  rate  enhanced 
by  the  presence  of  previously  cooled  Beryllium  ions. 

Then  beam  pulse  of  Xe'*'*^  ions  from  the  EBIT  (analyzed  by  a similar  magnet, 
passing  through  ion  optics  and  a bender  now  with  potentials  applied)  is  then  similarly 
captured  into  the  Penning  ion  trap  previously  filled  with  cold  Be”^  ions.  Electrode 
switching  and  collisions  with  the  xenon  ions  initially  heat  the  Be^  ions.  This  heating 
is  associated  with  an  increase  in  fluorescence  rate.  This  laser  scattering  quickly  re- 
cools the  Be^  ions,  which  cools  the  ~ 400  caught  Xe  ions.  The  mixed  ion  cloud  cools, 
different  charge-to-mass  ratio  ions  separate  with  the  cold  highly  charged  Xe  ions  in 
the  center  of  the  trap  and  the  Be^  ions  in  an  annulus  that  we  image  using  scattered  313 
nm  light.  A more  detailed  account  of  the  procedure  and  results  can  be  found  in  ref. 
[7].  Laser  induced  fluorescence  from  a magnetic  dipole  transition  in  the  ground  term 
of  a suitable  highly-charged  ion  can  be  potentially  used  to  image  the  highly  charged 
species  in  the  trap[4]. 
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FIGURE  1.  Setup  to  inject  MeVVA  ions  into  trap  at  LLNL.  Magnet  is  used  to  separate  charge  states 
produced  by  MeVVA. 

The  shape  of  Be^  fluorescent  image  can  be  used  to  determine  the  plasma 
parameters  such  as  density  and  rotation  rate.  An  ion  species  of  interest  can  be  excited 
by  cyclotron  resonance  excitation  (in  the  present  trap  configuration  the  frequencies 
associated  with  0)+,  the  modified  cyclotron  frequency  is  excited)  resulting  in  heat 
transfer  to  the  Be"^  ions  and,  typically,  an  increase  in  scatter  of  cooling  laser  light  by 
the  hotter  Be^  ions.  Both  the  frequency  and  time  response  of  this  heating  signal  are  of 
interest. 

A preliminary  study  of  cooling  laser  detuning  (cooling  power)  with  various  radio 
frequency  excitations  (heating)  of  an  ion  cloud  resulting  from  a Be^  - Xe^"*^  capture 
shows  interesting  features  (see  figure  3).  These  signals  may  possibly  indicate 
crystallization,  some  other  configuration/coupling  change  between  species  or  some 
systematic  effect  that  is  not  understood.  The  laser  tuning  was  monitored  (using  the 
fundamental  of  the  doubled  laser)  against  an  I2  absoqMion  cell  for  stability  in 
frequency  and  amplitude.  Similar  signals  have  been  assigned  for  small  numbers  of 
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ions  in  a Paul  trap[9].  Unfortunately,  pump  laser  failure  and  other  considerations 
prevented  further  studies.  We  wish  to  pursue  similar  studies  at  the  new  location  of  the 
experiment. 

The  trap  contents,  when  warmer,  can  be  examined  with  the  tuned  circuit  used  for 
electronic  cooling  as  described  above.  The  trap  contents  can  also  be  destructively 
probed  by  releasing  the  ions  into  various  charged  particle  detectors. 


FIGURE  2.  Image  of  Be^  ion  cloud  containing  Xe’'**and  daughter  charge  exchange  ions.  B field  is 
roughly  vertical.  With  the  cooling  laser  detuning  and  simulations  it  is  estimated  that  the  Xe  ion 
temperature  is  approximately  IK.  The  image  has  been  analyzed  to  give  an  Be^  ion  density  2.16  lO''*  m’^ 
and  a Xe’^^  density  of  6.4  lO’^  m'^.  This  corresponds  to  a F value  of  about  577. 
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FIGURE  3.  Cooling  laser  tuning  curves.  Main  perpendicular  cooling  laser  is  used.  P|  is  base  curve 
with  no  extra  radio  frequency  heating  applied.  F-rom  other  experiments  the  plasma  temperature  is  about 
1 Kelvin  when  the  the  detuning  for  the  PI  curve  is  200  Mllz.  The  other  curves  have  small  amounts  of 
radio  frequency  power  applied  near  the  Be*  cyclotron  frequency.  The  trap  contains  Be*  , Xe’’*  and 
charge  exchange  daughter  products.  Plasma  conditions  arc  similar  to  that  of  Figure  2.  The  same  ion 
cloud  is  used  for  all  curves. 


SIMULATION 

A molecular  dynamics  code  originally  developed  by  E.  L.  Pollock,  and  used  by  H. 
DeWitt  et  al.  to  model  strongly  coupled  plasmas  [10]  has  been  modified  to  run  on 
personal  computers.  It  was  modified  to  exclusively  calculate  the  interaction  of  two  ion 
species  confined  in  a Penning  trap.  Simple  heating  or  cooling  of  each  individual 
species  was  provided. 
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FIGURE  4.  A “snapshot”  x-y  projection  of  the  cloud  of  64  Be*  and  64  Xe'*'’*  ions  used  for  heating  rate 
calculations.  This  projection  is  for  0.5  Kelvin  axial  cooling  of  the  Be*  ions  and  1000  Kelvin/(particle  • 
second)  axial  heating  of  the  Xe''"'*  ions.  The  trap  parameter  is  q = -2.0  corresponding  to  about  540  Volts 
applied  to  our  trap.  The  rotation  frequency  of  the  cloud  is  approximately  865  kHz. 


Simulations  of  the  motion  of  5 to  64  multiply  charged  ions  (experimental  -100- 
500)  and  64  to  200  Be^  (experimental  -10^  ) ions  were  performed  using  fields, 
potentials,  ion  masses,  and  ion  charges  correspond  to  those  used  in  the  experimental 
work.  Initially,  the  ions  were  arranged  in  a body-centered  cubic  lattice,  with  random 
assignment  of  initial  locations  of  each  species.  The  center  of  mass  was  placed  at  the 
trap  center  to  simplify  calculation.  The  ions  were  assigned  pseudo-random  initial 
velocities  corresponding  to  an  initial  temperature,  with  the  net  linear  momentum  of  the 
ensemble  set  to  zero.  The  initial  ion  density,  aspect  ratio  of  the  lattice,  and  ion 
rotation  velocity  were  assigned. 
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FIGURE  5.  A “snapshol"  cylindrical  R-z  projection  of  the  ion  cloud  described  in  Figured. 

The  electrostatic  force  between  each  pair  of  particles  was  calculated  using 
Coulomb’s  Law  between  point  charges.  The  electrostatic  trap  potential  was  modeled 
as  an  ideal  hyperbolic  potential,  with  the  voltage  applied  between  ring  and  end  caps  of 
a trap  of  specified  size  characterized  by  a single  parameter  q,rap.  For  the  9.34  mm 
scale  length  trap,  q,„p  = -2  corresponds  to  about  540  V.  This  approximates  the  voltage 
at  which  the  axial  oscillation  of  a Be'^  ion  is  resonant  with  the  2.5  MHz  tuned  circuit. 
The  value  q^p  = -0.039  corresponds  to  a shallow  axial  well,  used  experimentally  for 
lower  ion  plasma  rotation  frequencies,  resulting  in  greater  interaction  stability  between 
the  laser  beam  and  the  plasma.  Dynamic  features  associated  with  ion  confinement 
measurements,  such  as  r.f.  excitations,  thermal  bath  coupling  to  tuned  circuits,  and 
image  charges,  were  not  modeled. 

The  equation  of  motion  of  the  ions  was  integrated  in  the  lab  frame  using  the  Verlet 
method  [11],  based  on  the  Stoermer  formula  rj(i  + t)  = 2Pj(t)  - ri(t  - x)  + x"ai(t)  + 0(x‘‘). 
Here  aj  is  the  acceleration  of  ion  i due  to  the  electric  and  magnetic  forces  of  the  trap. 
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as  well  as  the  Coulomb  forces.  Since  the  trap  symmetry  axis  is  aligned  with  the 
uniform  magnetic  field  B,  this  equation  is  simplest  to  apply  in  the  z coordinate.  For 
the  X and  y radial  coordinate,  the  time-centered  velocity  v-,(t)  = (ri(t  + t)  -rj(t  - t))/2t  -i- 
O(t^)  was  used  to  evaluate  the  Lorentz  force  in  the  Stoermer  formula  and  to  get 
implicit  expressions  for  x(t  + x)  and  y(t  -t-  x)  in  terms  of  present  acceleration,  x and  y 
position,  and  a single  past  time  step  of  stored  past  x and  y positions.  A more 
sophisticated  approach  is  described  in  the  work  of  Spreiter  and  Walter,  which  we  hope 
to  incorporate  in  future  simulations[12].  The  parameter  x was  chosen  to  be  a fraction 
of  the  period  of  the  initial  plasma  oscillation,  so  that  it  is  also  small  enough  to  follow 
the  largest  cyclotron  frequency  of  the  simulation.  It  was  found  that  small  values  of  x 
(~30  - 300  ps  compared  to  -150  ns  Be^  cyclotron  period)  resulted  in  a stable 
simulation  start. 

The  mean  squared  axial  velocity  component  (proportional  to  axial  temperature  Tj 
of  each  species)  is  averaged  for  every  other  step.  Every  1000  steps,  an  estimate  of  the 
rigid  body  rotation  frequency  was  made,  and  in  that  reference  frame,  the  mean  squared 
radial  velocity  for  each  species  (proportional  to  T,)  was  calculated  and  recorded,  along 
with  the  fixed  frame  kinetic  energies  and  potential  energies  of  the  ions.  When 
simulating  cooling,  if  the  running  average  axial  temperature  of  an  ion  species  exceeds 
the  specified  temperature  by  3 %,  the  former  z-axis  positions  of  the  species  being 
“cooled”  (proportional  to  axial  velocity  in  a harmonic  well)  were  re-scaled  using  a 
Taylor  series  expansion  to  approach  the  desired  temperature.  The  other  ion  species 
could  also  be  “heated”  axially  by  a similar  scaling,  allowing  a fixed  amount  of  kinetic 
energy  to  be  added  on  average  to  each  ion.  The  temperature  of  the  heated  species  can 
approach  a steady  state  with  the  cooled  species.  Typically,  the  target  temperature  is 
lowered  in  steps  as  the  system  approaches  quasi-steady  state. 

Major  results  of  the  simulation  include  observation  of  the  onset  of  properties  of  a 
highly  magnetized  non-neutral  plasma  [13,14].  The  multiply  charged  ions,  with  higher 
cyclotron  frequencies,  exhibited  this  effect  at  higher  temperatures  than  the  Be'^  ions, 
(ranging  from  40  K to  2 K for  Xe'*'*^)  depending  on  density  (or  equivalently  cloud 
rotation  frequency),  as  the  “radial  cyclotron  temperature”  decoupled  at  an  exponential 
rate  from  the  “axial  temperature”. 

For  example,  a simulation  of  64  Be^  ions  and  64  Xe  ion  in  a qirap  = -2.0  well  and 
with  a cold  rigid  body  rotation  frequency  of  865  kHz,  this  temperature  was  about  0.5 
Kelvin.  A snapshot  of  this  simulation  is  shown  in  Figures  4 and  5.  The  ions  are 
separated,  as  discussed  below.  However,  the  radial  and  axial  motions  remain  coupled 
through  the  ion  cloud  rotation.  The  Be^  ions  were  not  significantly  magnetized  for 
simulated  temperatures  > 1 Kelvin. 

Typically,  the  multiply  charged  ions  and  Be^  ions  centrifugally  separate,  as  Figures 
3 and  4 show.  This  centrifugal  separation  occurs  when  the  plasma  rotates  as  a rigid 
body.  The  “pseudopotentials”  associated  with  the  plasma  rotation  for  each  species 
then  differ,  because  of  the  different  mass-to-charge  ratios,  causing  the  ions  with  higher 
mass  to  charge  ratio  (in  this  case  Be^)  to  move  radially  outward  and  separate  from  the 
HCIs.  The  separation  temperature  was  found  to  depend  on  the  trap  parameter  qtmp 
and  rotation  frequency  (bounds  of  which  depend  on  q,r„p)  as  expected,  although  the 
current  code  (no  spacial  averaging  / small  numbers  of  particles)  and  mode  of 
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simulation  (step  down  of  set  tempefature)  does  not  track  the  separation  in  detail. 
Simultaneously  ordering  of  the  HCIs  can  be  seen,  first  as  shells,  then,  at  lower 
temperatures,  rings  of  ions  are  apparent.  The  shell  surfaces  then  have  the 
characteristic  triangular  structure  of  2D  Coulomb  systems  [14,15]. 

Additional  work  involves  the  study  of  the  effect  of  heating  sources  on  the  multiply 
charged  ion  temperatures  given  an  experimentally  measurable  Be"^  temperature.  From 
simulated  motion  and  atomic  data,  spectroscopic  lineshapes  can  be  calculated  and 
parameters  for  future  experiments  estimated.  For  example,  Ca''*^  has  been  caught  and 
electronically  cooled  in  our  trap.  Unfortunately,  Be^  laser  cooling  was  not  functioning 
at  that  time  due  to  laser  failure.  Calculations  are  proceeding  to  model  this  system  to 
optimize  our  procedure  for  a future  experimental  opportunity. 


THE  MOVE 

The  Retrap  facility  has  been  moved  to  the  Lawrence  Berkeley  National  Laboratory. 
The  trap  is  now  sited  on  the  roof  of  a target  cave  at  the  LBL  88”  cyclotron[16].  The 
system  is  now  undergoing  preliminary  testing  with  a MeVVA  ion  source.  Several 
electron  cyclotron  resonance  (ECR)  ion  sources  will  be  available  to  continue  work 
with  highly  charged  ions.  The  laser  capabilities  needed  for  cooling  and  spectroscopy 
is  to  be  reestablished  in  the  near  future.  The  cyclotron  and  associated  ion  source 
technology  will  allow  the  trapping  and  study  of  rare  isotopes.  The  effect  of  charge 
state,  for  example,  on  various  nuclear  decay  mechanisms  can  be  studied. 
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